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Abstract 

We discuss the QCD constraints on e + e~ — > baryon-anti-baryon 
close to threshold, in light of the puzzling experimental data which 
indicate that close to threshold a(e + e~ — ► nn) > a(e + e~ — > pp). We 
focus on the process e + e~ — ► AA, which is particularly simple from 
the theoretical point of view. In this case it is possible to make exact 
QCD predictions for the relative yields of the four members of the A 
multiplet, modulo one crucial dynamical assumption. 
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It has been pointed out recently that the experimental results 0-|| 
indicating that 

cr(e + e~ — > nn) > a(e + e~ — > pp) (1) 

over the threshold region, W = 1.88 4- 2.45 GeV, are rather surprising. 

Indeed, if the fragmentation amplitude F(u — > proton) = F(d — > neutron) 
is not smaller in magnitude than F(d —>■ proton) = F(u — > neutron), then 
the pp cross section is larger than that for nn. This is so since the amplitude 
for uu pair production by the virtual photon is enhanced by the ratio of 
quark charges \Q U /Qd\ — 2, relative to the corresponding dd amplitude. 

It has therefore been argued that such reasoning - inspired by 

perturbative QCD, does not apply at NN threshold. The q — > B and q — > 
B fragmentation functions are not independent, annihilation channels into 
multipion states dominate, and other, say Skyrme-like descriptions for the 
nucleons are more appropriate than the "naive" three quark picture. 

Let us briefly restate the main argument of @. The 1 = 1 (say p,p',p", 
etc.) intermediate states leading to an even number of pions dominate the 
annihilation into multi-pion states. This is expected from the direct quark 
EM couplings of the uu and dd quarks to the photon in an ideal "nonet" (or 
"quartet" when ss production is neglected) symmetry of the vector multiplets 
(we would like to emphasize that such symmetries arise when purely gluonic 
1 1 = channels are neglected). 

If further the 1 = 1 amplitudes dominate the NN 1 annihilation chan- 
nel - as the ratio of pp — ► K + K~ / K°K° suggests, then 1 = 1 intermediaries 
dominate NN —>■ e + e~ at threshold. 

If \A(I = 1)| > \A(I = 0)| then a(e + e~ -> pp) = a(e + e~ -> nn). Further 
admixture of some 1 = amplitude can yield the desired experimentally mea- 
sured cross section ratio for a range of the parameters e = \A(I=0)/A(I=1)\ 
and ArgL4(I=0)/A(i=l)]. 

In the following we focus on another baryon-antibaryon threshold process 
which was briefly alluded to in Jl|, namely the (pair) production of / = 3/2, 
J = 3/2 A states: A++, A+, A and A - . 

Our discussion does not contradict the above phenomenological consider- 
ations of the NN threshold process for the 1 = 1/2 nucleons. Yet it suggests 
a complementary point of view and possible "exact" QCD predictions even 
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in this non-perturbative regime, modulo one crucial assumption, as discussed 
in detail below. 

We start by comparing the two processes, 

e + e~ — ► A ++ A++ : e + e~ — > 7* — > uu — > mm + uuu (2) 

and 

e + e- -> A" A - : e+e" -> 7* -> dd ^ ddd + ddd (3) 

a- 

Unlike the e + e~ — ► iViV case, where one has a mixed-symmetry nucleon 
with two different fragmentation functions, only one quark-baryon fragmen- 
tation function is involved here, namely F(u — > A ++ ) = F(d — > A - ) . Thus 
we have a clear prediction of the "naive" quark model, 

a(e + e- -> A++A++) = Aa(e + e~ -> A^A - ) (4) 

The key point that we wish to emphasize is that this prediction is far 
from naive. Eq. ([|) above and other predicted ratios of reactions involving 
A(1238) baryons become exact when the following assumptions are made: 
(a) we take the m u = m d limit; (b) we work to lowest order in c^qed ] ( c ) we 
consider only "flavor connected" processes in which the initial q^i produced 
via the initial (virtual) photon do not annihilate into gluons and appear in 
the final hadrons. 

While (a) and (b) are standard, assumption (c) is not as clear-cut. It 
entails a l/N c approximation extended to baryons ||-[]nj, and in the 
present context it is equivalent to the famous Zweig (OZI) rule. [] 

The 1/N C factor suppresses flavor-disconnected contributions where the 
initial external current couples via an internal quark loop, relative to flavor- 
connected processes where the initial quarks produced by the EM current 
are incorporated into final hadrons. 

We wish to stress that the issue of OZI rule applied to baryons is quite 
delicate: while it works well for mesons, it is well known that in the baryon 

^The extension to baryons of involves some subtleties which will be discussed in some 
detail a bit later. 
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sector there are large OZI violations, at least in one case, that of ss meson 
production in BB annihilation at rest fll"l, ||i~2]| . It is important to see whether 
this is compatible with assumption (c), to be discussed in more detail later. 
Consequently, we now proceed on the basis of all three assumptions. Clearly, 
confronting the ensuing predictions with experiment will be a sensitive test 
of assumption (c). 

First, we note that relation (|J) is valid to all orders in perturbation theory. 
For any QCD Feynman diagram contributing to e + e~ — > A ++ A++, there is a 
corresponding diagram, with all w-quark propagators replaced by (identical!) 
d-quark propagators, which contributes to e + e~ — > A~A~. Because of the 
flavor independence of gluonic couplings these two diagrams will therefore be 
identical as far as QCD is concerned. The only difference is that Qd = —1/3 
replaces Q u = 2/3 in the coupling to the external photon. Thus we expect 
that the amplitudes for the two different processes will satisfy 



A(e + e- -> A++A++) = -2A(e+ e - -> A" A - ) (5) 

We omitted all the kinematical variables, namely the external momenta 
and helicities, which should be the same for the two processes. 

Upon adding all the squares of the helicity amplitudes we obtain the 
differential cross section and above follows. 

Needless to say, we do not use here any approximate (S-wave) symmet- 
ric, naive, quark model or any other model for the wave- (or fragmentation) 
functions, and allow for any admixture of orbital/radial excitations - and 
contribution of states with any number of gluons and/or extra qq pairs. The 
only issue at stake here is the u «-> d exchange symmetry. This is a vectorial 
flavor symmetry. The Vafa-Witten theorem ]nj, which relies on rigorous 
QCD inequalities (see |l4j and [15|] for reviews) states that such symme- 
tries do not break down spontaneously, and hence the all-order perturbative 
relations should be true in the full theory. 

Our argument is actually independent of any perturbative expansion. All 
hadronic amplitudes can be obtained by Fourier transformation and ana- 
lytic continuations of Euclidean n-point correlators of (color singlet) local 
observables. 

The relevant correlators for the processes of interest are the three-point 
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functions of the form: 

K A++ (x,y,z) = (0\J EM (x)A ++ (y)A++(z)\0) (6) 

with J^em(x) = |ti(x)7 At u(x) — |d(x)7 /1 d(x) + . . . the electromagnetic cur- 
rent and A ++ (x) a local operator trilinear in quark fields with the quantum 
numbers of the A ++ particle, 

A ++ (x) = u a ,{x) u bJ (x) u c , k (x) e abc V> k (7) 

with a,b,c (i,j,k) color (spinor) indices, e and T guarantee coupling to an 
overall color-singlet baryon with total spin J = 3/2. We have an analogous 
expression for Ka-(x, y, z). 

There are two contraction patterns of the fermionic operators in the above 
correlators - the "flavor connected" and the "flavor disconnected" . 

In the first case the initial quark ( ant i- quark) created from the vacuum 
via Jem at x propagates to the final y (or z) vertex. In the second case, 
these quarks annihilate, forming a closed loop which starts and terminates 
at the same point x. The path integrals for evaluating these contributions 
are respectively 

K c ^ ected (x, y, z) =Q U J d f i(A)Tr{^(x,y)eTS^z,y)S^z,y)e^S^z,x)} 

(8) 

and 

K d £T n ix,y,z) = 

£ Q t J d^A)Tr{^Sf(x,x)} Tr{ers£(z,y)S£(z,y)S£(z,y)eI*} (9) 

i=u,d,s 

where the above traces are in color/spinor space , Sf(x,y) is the propagator 
from x to y of the i-th flavor quark in a given background gauge field A^x), 
and 

dfi(A) =VA^x)e-S d4 < E2+B2 ) UDetipA + rm) (10) 

i 

is the measure in the path integral. 

The two contraction patterns are illustrated schematically in Fig. 1(a) 
and 1(b), respectively. 
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eT y eTz eT y 




(a) (b) 

Figure 1: The two contraction patterns for the fermionic operators determining 
the three-point point functions (JemAA) in eq. ftfy: (a) - the flavor- connected 
contraction, eq. (b) - the flavor-disconnected contraction, eq. (|^). The quark 
loop in (b) includes a sum of contributions from all the light flavors. 

Baryon anti-baryon production is expected to be strongly suppressed in 
the large- N c limit. Still, it is meaningful to compare the N c hierarchy of the 
two classes of processes - 1(a) and 1(b). 

Apart from an external color-singlet insertion, 1(a) is essentially the same 
(in N c terms) as the baryonic propagator. The contribution 1(b), on the 
other hand, has a quark loop in addition to the baryonic propagator. As 
is well known from the double (single) line counting rules for gluonic(quark) 
propagators JFj, such a quark loop introduces an extra g 2 factor, yet it does 
not increase the number of color traces. 
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Hence the 1(b) flavor-disconnected contraction is suppressed by 1/N C rel- 
ative to the connected contraction 1(a). Analogous expressions can be writ- 
ten for the correlator K&-(x,y,z) and for the flavor-connected and flavor 
disconnected-parts thereof. 

This argument is quite simple for mesons, but the extension to baryons, 
although yielding the same qualitative conclusion [pf-|]T0|, involves some sub- 
tleties. 

Although the result is not new, for the reader's convenience we review the 
argument in the following. Consider the large-iV c baryon propagator with 
the leading quark-loop contribution. It is the diagram with a quark-loop 
correction to a gluon propagating between two valence quarks. There are 
two vertices where the gluon couples to the valence quarks and two vertices 
where the quark loop couples to the gluon. Each vertex is associated with 
1/y/Nc, and thus each individual diagram of this kind is suppressed by 1/N% 
with respect to the diagram with only valence quarks propagating. 

The subtlety arises because there are N C (N C — l)/2 such pairs of valence 
quarks in the baryon, and so the total contribution of the diagrams with a 
quark loop is of order 0(1) and thus appears to be unsuppressed. 

However, as stressed by Witten ||, diagrammatic expansion is not a 
very convenient way to study the large- N c limit of baryons, because in the 
perturbation series each successive order grows as a higher power of N c . Thus 
the diagram with one gluon exchanged between the two quarks in the baryon 
scales like 0(N C ), (two vertices ~ l/yfN~ c each, multiplied by combinatorial 
factor of N C (N C — l)/2 of possible quark pairs), a diagram with two gluons 
exchanged between two different pairs of quarks is of order N%, etc. This 
divergence is simply an an artifact of the baryon mass being ~ 0(N C ). 

The crucial point for our discussion is that the quark loop is 1/N C sup- 
pressed with respect to the contribution from the same diagram without the 
quark loop: thus the diagram with one gluon exchanged between two valence 
quarks scales like 0(N C ), while the same diagram with an additional quark 
loop correction to that gluon scales like 0(1), etc. 

The general arguments given by Witten make it clear that the total con- 
tribution of gluons to the baryon propagator is of the same order as that of 
the valence quarks, i.e. linear in N c . Since any diagram with a quark loop is 
l/N c suppressed with respect to the same diagram with only valence quarks 
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and gluons, it follows that quark loops are suppressed by l/N c with respect 
to the valence quarks P-fL0| . 



In the real world, this suppression manifests itself, in particular via the 
absence of " exotic" resonances (involving extra qq pairs) and the narrowness 
of the ordinary, non-exotic resonances. These hold equally well both for 
mesons and for baryons. 

Apart from the overall quark charge factors of | and -~, respectively, the 
connected K&++(x,y, z) and K&-(x,y, z) are identical. This identity holds 
in a very strict sense, namely pointwise in the path integrals and @, i.e 
for each gauge field configuration A^ix) separately. Hence it is guaranteed 
to hold also for the integrated quantities. This is true in any scheme, such as 
lattice gauge theory, used in order to define and regularize the path integral 
in a rigorous way. 

Thus we conclude that 

K%?r cted (x, y, z) = -2K^ nected (x, y, z) (11) 

holds for any x, y and z. It will therefore hold for the physical amplitudes for 
e + e~ annihilations into the corresponding AA baryons which are obtained by 
common manipulations of the two Euclidean correlators. Thus, to the extent 
that we neglect the flavor-disconnected part, relations (f|) and @ ensue.0 

It is important to note that for the case at hand the X conn ' contribution 
is suppressed not only by the 1/N C . There is also an independent effect which 
works in the same direction, namely a partial cancellation of the contributions 
of the three light quark flavors u, s and d. The cancellation, due to Q u + 
Qd + Qs = 0, is exact in the SU(3) flavor symmetry limit. The physical 
world is quite well approximated by this limit, since the m s — m u /d mass 
difference is negligible in comparison with the mass scale relevant for the 
problem ~ = 1 GeV. In this case the three terms, due to u, d and s 
loops appearing in K dlsconn ' above are identical as far as QCD is concerned.^ 

* These predictions are shared by another frequently used nonperturbative approach. 
This is the model utilizing qq pair production via tunneling in a chromoelectric flux tube 
[fl6| , incorporated in the Lund model (l7|] . Flavor-disconnected parts are neglected in such 
models. 

*It should be emphasized that unlike the "universal" 1/N C suppression, this charge- 
cancellation mechanism is relevant only in the specific energy range considered here. At 
much lower energies only uu and dd contribute, and at much higher energies cc contribute 
as well. 
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The above arguments apply also for the case of annihilations into NN 
states. However, in the latter case we have, even when only the flavor- 
connected parts are retained, two distinct contributions. One contribution 
occurs when the primary quark produced by the photon is one of the two 
same-flavour quarks in the nucleon (namely the u in the proton or the d in 
the of the neutron). The other contribution corresponds to the case when the 
other quark (d for proton or u for the neutron) is the primary one. There is 
no a priori model-independent relation between these two contributions, nor 
do we know of any systematic expansion like 1/N C which could help relate 
them, and hence no analogous firm prediction can be made here. 

In passing we note that also 

KflT n ix, y, z) = Kt s - conn ix, y, z) (12) 

and the resulting flavor-disconnected contributions to the physical ampli- 
tudes are also identical. This is indeed expected, as the disconnected am- 
plitudes correspond to purely gluonic and hence 1 = intermediate states, 
which couple equally to all members of the A(1238) multiplet. 

Considering next the 1=1 and 1 = isospin channels we find that 
eq. ( |TTD above implies that for the 7* — > A A, 

A(I = 1) = 2A(I = 0). (13) 

Interestingly, this is consistent with the dominance of the 1=1 channel 
in e + e~ — > NN, as inferred from phenomenological analysis of the data in 
Ref. [pj. It is reassuring that here it follows directly from QCD, augmented 
by assumption c. 

Isospin symmetry and the above amplitude ratio imply 

A(A ++ ) : A(A+) : A(A°) : A(A~) = 2 : 1 : : -1 (14) 

Thus the e + e~ — > AA cross section ratios 4:1:0:1 discussed in ffl, are in 
fact exact QCD results, modulo assumption (c) above, i.e. neglect of the 
"flavor disconnected" parts.[] 

^For the processes e + e~ — > NA, involving nucleon and delta production, only the 1=1 
amplitude contributes and the ratios are fixed by isospin alone. 
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It is interesting to consider the effect of possible resonances in the 1 = 1 
(p, p', p") and the 1 = channel (u, uj', uj"). Clearly, if at any particu- 
lar energy an I = or an / = 1 resonance dominates, the above ratio of 
A(I = 1)/A(I = 0) cannot be maintained. 

We note however that precisely in the large N c limit, when the quark 
annihilation diagrams are neglected, the 1 = 1 and 1 = states are degener- 
ate: m w = m p , m w i = m p >, etc. In reality, it is much easier to maintain the 
A(I = 1)/A(I = 0) ratio, since the resonances in the W > 2 GeV region are 
broad and overlapping. 

Independently of the ordinary p^ and U)( n > resonances and/or annihila- 
tion channels, dynamical enhancement of the gluonic 1 = intermediate state 
could result if a relatively narrow glueball state with a vector 1 quantum 
numbers occurs somewhat above the A(1238)A(1238) threshold. It would 
manifest itself via a transient deviation from (|^) when the resonant energy 
is traversed. 



Qualitative arguments JT3| suggest that the lightest 1 state in the 
"three-gluon" sector is at least 50 % heavier than the lightest scalar glue- 
ball in the "two-gluon" sector, with mo++ — 1.6 4- 1.7 GeV. In other 



words, mi — <S 2m a- Lattice calculations |y| indicate a much larger value, 
m 1 — ~ 3.8 GeV. If this is indeed the case, the 1 glueball is too far above 
the AA threshold. It would be very interesting to settle this issue. 

Following the suggestion of Refs. |IJ and 0, we next consider the 77 — > 
BB processes. We again compare A ++ and A~ pair productions. If we 
neglect the disconnected part, then the same arguments imply a ratio of 
4 of the flavor-connected amplitudes and 16(!) for the corresponding cross 
sections, as discussed briefly in 0. 

It should be emphasized, however, that the neglect of the flavor- 
disconnected part relies here, unlike for the previous single-photon case, on 
the 1/N C suppression only and not on cancellation of charges as well. There 
is no such cancellation here, as only the squares of the charges enter. 

In addition, in the two photon case we have also the 0~ rj anomalous 
channel, which implies that the flavor-disconnected parts may be enhanced. 
However, as we move even slightly above the BB threshold, the number of 
partial waves contributing increases rapidly and the relative importance of 
the pseudoscalar channel decreases. 
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The present work is indirectly triggered by present data on e + e~ — > nn/pp 
at threshold, which seem to confound the naive QCD expectations. Our 
main point is to note that there is a similar process, in which the theoretical 
analysis is much more straightforward, namely e + e _ — > A A. In this case 
QCD, augmented by a well-motivated dynamical assumption, makes clearcut 
and largely model-independent predictions. Future experimental tests of 
these predictions will be a sensitive probe of this dynamical assumption. 
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